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ABSTRACT 

We calculate the evolution of a low-mass (M < lO^M©) spherically symmetric 
density perturbation in the Oft/i^ = 0.02, Qm = 0.35, = 0.65, h = 0.72 Universe. 
The results are compared with the ones that assume no cosmological constant and 
the flat, dark matter dominated Universe. We include thermal processes and non- 
equilibrium chemical evolution of the collapsing gas. We find that direct formation of 
bound objects with such masses by z = 8 is unlikely so in fact they may form only 
through fragmentation of greater objects. This is in stark contrast to the = 1 pure 
CDM cosmology, where low-mass objects form abundantly at redshits z > 10. 

Key words: hydrodynamics - instabilities - dark matter ~ early Universe. 



1 INTRODUCTION 

Present observations of the C MB an isotropies by the 

|2002[ ), MAX IMA- 1 Cee 



et al 



BOOMER ANG ( iNetterfield et al 
|200l| ), CBI ( JPadin et al 



200l|)and DASI (Prykeetal 



2002) experim ents combin e d wit h the data from the Supe r 
(|199^ 



nova project (Riess et al 



Perlmutter et al 



(1999)) 



suggest that although the geometry of the Universe is flat, 
its matter content is dominated by the vacuum energy or so- 
called quintessence. The background cosmology has a very 
strong influence on the structure formation. In this paper 
we address the formation of first bound objects in the flat 
ACDM cosmology. 

The structure formation on both large and small scales 
has been studied very carefully in the dark matter domi- 
nated cosmology, in particular in the flat pure CDM model. 
This is not so for the vacuum-energy-dominated models, 
which are of prime relevance according to the above obser- 
vations. Perhaps one of the most interesting objects to study 
are the first bound objects in the Universe - the first quasars 
and the Population III stars that were the first sources of 
photons after the light of CMB faded away, and ended the 
'Dark Ages' of the Univ erse. Some of thes e calcu lations were 
described in a review by Barkana fc Loet ( 2001 ). The recent 
identific ation of quasar s at z ~ 6 in the Sloan Digital Sky 
Survey (Fan et al. 2001) strongly motivates such an investi- 
gation. 

Our code is similar to the one by Thou l & Weinberg 



( tl995| ) and by [Haiman, Thoul fc Loebj ( |1996D . We consider 
the evolution of a single spherically symmetric density per- 
turbation in the early Universe starting soon after recombi- 
nation until it finally forms a bound stationary object. Our 



aim here is to extend their analysis, which was restricted to 
the = 1 CDM cosmology, to the A-dominated Universe 
and compare the results with those for the fiat, dark matter 
dominated Universe but with more up-to-date cosmological 
parameters [h = 0.72, Ot/i^ = 0.02). 

We start tracing the initial expansion of the perturba- 
tion at high redshift when its density contrast is still in the 
linear regime. Then we follow decoupling of the perturba- 
tion from the Hubble flow and its subsequent collapse and 
formation of a virialized cloud. We include the gas dynamics 
and various cooling and heating processes operating in the 
expanding and collapsing cloud. The chemical evolution of 
the collapsing primordial gas cloud is also accounted for. 

After the initial collapse, a virialized gas cloud is 
formed. The kinetic energy of the infalling gas is dissipated 
through shocks and the cloud becomes pressure-supported. 
We study the further evolution of the cloud which is deter- 
mined by its ability to cool sufficiently fast. The most impor- 
tant cooling mechanism for low-mass clouds is the radiation 
of excited H2 molecules. The presence of a small amount of 
the molecular hydrogen H2 is crucial for triggering the final 
collapse of such clouds which could form the first luminous 
object in the Universe. 



2 NONLINEAR EVOLUTION OF A 

SPHERICALLY SYMMETRIC DENSITY 
PERTURBATION 

A major approximation we use is the assumption of spherical 
symmetry. This assumption is best justified in the case of 
first objects formed in the Universe which are supposed to 
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originate from the rare highest fluctu ations in the primordial 



density field (Bardeen et al. 19861 ). It allows us to focus 
mainly on gas dynamical processes which are expected to 
control the collapse of low-mass clouds which are supposed 
to form the first luminous objects. 

The spectrum of density fiuctuations has more power 
on small scales, hence the first nonlinear structures are ex - 
pected to occur on relatively small scales ( Vliller et al, 1999). 
As we are interested in small scales, much lower than the 
horizon, we use Newtonian gravity and treat the expansion 
of the Universe as a hydrodynamical flow. To describe the 
evolution of a spherically symmetric density perturbation in 
the nonlinear regime we use Lagrangian coordinates. We di- 
vide both baryonic and dark matter into concentric shells. 
The dynamics of primordial matter is given by the equations 
below. 

The continuity equation for baryonic matter reads 

— — =4TTrtgb, (1) 
an 

where rt is the radius of a sphere of mass Mi,. The radial 
velocity of the surface of this sphere is 

drb 

and the acceleration satisfles the dynamical equation 

dvb _ _ 2 dp ^ GM{rt) 
It ^ ~ "IM: 7i ' 



(2) 



(3) 



where M{rb) = Mb{n) + Mi^{n) + 2{M,.ad{rt) - MAin)) 
is the total mass within radius rt including the effects of ra- 
diation and vacuum energy components of the stress-energy 
tensor. Two last contributions result from the CMB and cos- 
mological constant components of the energy density of the 
Universe, 



Mradirt) 

and 

MAin) = 



Hi 
2G 



i^rad(z + 1) rt 



m 

2G 



(4) 



(5) 



where Qrad = SnGaTc m b''' / {3Hq c^) is the present contribu- 
tion of the CMB to the total energy density of the Universe 
in t erms of the critical density . Alth ough, as we have shown 
in ( ^tachniewicz fc Kutschera 2001), at least the radiation 
component is not negligible, both dark energy and radiation 
terms are not very important. 

The energy conservation condition for baryonic matter 
reads 

du _ p_dgb ^ ^cooi 
dt p'i dt 



Qb 



(6) 



where u is the internal energy per unit mass, p is the pressure 
and Qb is the baryon density. The last term in the eq.(^ 
describes cooling/heating of the gas, with Acoo! being the 
energy emission (absorption) rate per unit volume. 
We use the equation of state of the ideal gas 



f = (7 - ^)QbU, 



(7) 



where 7 = 5/3, as the primordial baryonic matter after re- 
combination is assumed to be composed of monoatomic hy- 
drogen and helium with a small (but very important) ad- 
mixture of molecular hydrogen H2. 
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Figure 1. The 'Sinusoidal' profile. 

Dynamics of dark matter is simpler, as we assume it to 
be collisionless. The continuity equation is 

dMirr, 



dra„ 



(8) 



where Vdm is the radius of dark matter sphere of mass Al^n 
The radial velocity of this sphere is 



dr, 



dm 



dt 



— Vdm, 



and the acceleration reads 

m ) 

~dr^ Til • 



(9) 



(10) 



To solve the above equations we must specify the cool- 
ing/heating function Acooi and the the initial conditions. 



3 INITIAL CONDITIONS 

Let us begin with the initial conditions. We start to follow 
the evolution of the perturbation soon after the recombina- 
tion at a sufficiently high redshift that the perturbation is 
still in the linear r egime. We choose th e initial redshift to 
be Zi = 500 like in 



Haiman et al. 



(1996) 



We apply the initial density profiles in the form of a 
s ingle spherical Fourier mode used also by Haiman ot al 
( [1996D 



Qi(r) = Q,iQc{l + 5i 



sin kr 

kr ' 



(11) 



where i = b, dm and Qc is the critical density of the Universe, 
Qc ~ 3H^ /8itG with H being the actual value of the Hubble 
parameter. The quantities St and 5dm measure, respectively, 
the baryon and dark matter density enhancement with re- 
spect to the mean densities Qb = ilbQc and Qdm = ^dmQc- 

For this profile there exist two distinguished values of 
the radius, _Ro and which correspond, respectively, to the 
first zero and the first minimum of the function sin(fcr)/fcr. 
They are shown in Fig. |^. 

Inside the sphere of radius _Ro ~ T^/k which contains 
mass Mo, the local density contrast is positive. The mass 
Mo and the radius _Ro will be referred to as the cloud mass 
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Figure 2. Comparison of the 'Sinusoidal' profile and the Gaus- 
sian one. 



and the cloud radius, respectively. The local density con- 
trast is negative for > r > Rq, with the average density 
contrast vanishing for the sphere of radius Rz = 4.49341/fc 
with the mass M^. According to the gravitational instabil- 
ity theory in the expanding Universe, the shell of radius 
Rz will expand together with the Hubble flow not sufi^ering 
any additional deceleration. This is why we regard this pro- 
file as very convenient in numerical simulations because it 
eliminates numerical edge effects and mentioned shell simply 
follows the Hubble expansion of the Universe. It can thus be 
regarded as the boundary of the perturbation and the mass 
Mz will be referred to as the bound mass. 

It is worth to note that for radii not greater than 3/4i?o 
this profile is very similar to the Gaussian profile 



Table 1. RMS overdensities for ACDM and flat CDM. 



gi{r) = D,tQc 



1 + Si exp 



2Rj 



(12) 



with Rf — l/2Ro. Both profiles are compared in Fig. |^. 

As the initial velocity we use the Hubble velocity for 
baryon matter. 



Vbir) = Hr, 



(13) 



w hereas for dark matter the initial velocity is (Haiman et 
' ]l99j ) 



al 



Vdmir) = Hr{l — g < 5dm >r). 



(14) 



The expression in brackets indicates averaging over the 
sphere of radius r. Slower expansion of dark matter at 
Zi — 500 results from the fact that the dark matter per- 
turbations start to grow earlier than do the baryon matter 
ones. 

Finally, the amplitudes of baryon and dark matter per- 
turbations, Sb and Sdm should be specified. To do so one 
should use the power spectrum corresponding to the initial 
redshift Zi = 500 that provides rms density fluctuations at a 
given mass Mq. For Cold Dark Matter model we use trans- 
fer functions calculat e d usin g the CMBFAST program by 
Seljak & Zaldarriagc (1996). The actual value of the rms 
density fluctuation is equal to 
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(15) 



where P{k) = 4n dl^^^lOOk"0. 05^ k^tf{k, zf . The value 
of normalization constant dnorm may be also obtained from 
the CMBFAST program. 

The rms fluctuation a{R,z) depends on the shape of 
the window over which the density is aversiged. However, 
this dependence is quite weak. In the table qI we list some 
values for z = 500 obtained by applying this formula to 
the transfer function and COBE-normalized normalization 
obtained from CMBFAST to the Gaussian profile of filter 
radius Rf — l/2Ro. Instead of Rf, we use Mo which is equal 
to the total baryonic mass within a sphere of the radius Rq. 

It is worth to note that these results are not much dif- 
ferent from t he res ults obtained by applying the fit by 



Eisen- 



stein & Hu ( ll99£| ). However, our results are systematically 



greater by about 23-24% for the A-dominated model and by 
about 7-10% for the pure CDM model. Probably it is due 
to the fact that the mentioned fit is less accurate for larger 
redshifts (especially for z > 30) and greater Qb/i^M ratios. 

As the baryon density perturbations start to grow only 
after the recombination, we have decided to set its initial 
value to Sb,i — 0.1Sdrn,i (Haiman et al. 1996). Let us note 
that CMBFAST predicts Sb,i to be about 0.2Sd,n,i in both 
A and da rk matter dominated model s. How ever, as we have 
shown in [Stachniewicz 



Kutschera 



(2001), if Sb,^/Sdm,^ is 



0.2 or less the results are pretty much insensitive to this 
ratio, runs with its values equal to 0.2, 0.1 and 0.0 were 
almost indistinguishible. 



4 BARYONIC MATTER COMPONENTS, 

CHEMICAL REACTIONS AND THERMAL 
EFFECTS 

To specify the cooling/heating function Acooi in Eq.(^ one 
should include all relevant thermal and chemical processes 
in the primordial gas. There are many papers discussing the 
mos t important contributions t o the function Acooi - see 
e.g. KatZj Weinberg, Hernquisj ( 1996 ). Al l formulas may be 
found in stachniewicz fc Kutschera| (2001). 

The primordial gas consists of neutral atoms and 
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molecules, ions and free electrons. In this paper we have 
taken into account nine species: H, H~, H^, He, He"*", He"*"^, 



H2, HJ and e". Some authors (e.g. GaUi fc Palla (1998)) in- 



clude also deuterium and lithium. 

The abundance of various species can, generally, change 
with time as the chemical reactions between species oc- 
cur and the ionization and dissociation photoprocesses take 
place in the hot gas. The chemical reactions include such 
processes as e.g. the ionization of hydrogen and helium 
by electrons, the recombination of ions with electrons, the 
formation of negative hydrogen ions, the formation of H2 
molecule s, etc. The full list of relevant che mical reactions is 
given in Stachniewicz & Kutschera (2001). Photoprocesses 
include ionization of neutral hydrogen H, helium He and 
He"*" by photons, and dissociation of negative hydrogen ions 
H~ and H2 molecules by photons. 

Time evolution of the number density of the component 
Hi is described by the kinetic equation; 



~dt 



9 



(16) 



The first component on the right-hand side of this equation 
describes the chemical reactions and the other one describes 
photoionization and photodissociation processes. The coef- 
ficients kim are reaction rates, the quantities pho- 
toionization or photodissociation rates and aimi and bji are 
numbers equal to 0, ±1 or ±2 depending on the reaction. 
All the reaction rates, photoionization and photqdisso cia- 
tion rates are given in Stachniewicz fc Kutschera ( 2001 ). 

The cooling (heating) function Acooi includes emissiv- 
ities (absorption rates) due to such processes as the col- 
lisional ionization of H, He and He^, recombination to 
H, He and He^, the coUisional excitation of H and He^, 
Bremsstrahlung, the Compton cooling and the cooling by 
deexcitation of H2 molecules. The formulae for the heat- 
ing/cooling contributions of vari ous processes are given in 
Stachniewicz & Kutschera (2001). 



NUMERICAL CODE AND INITIAL 
CONDITIONS 



The dynamical equations ([lD-(|lO[) are solved numerically. 
At each timestep also the chemical composition of the gas is 
updated by solving Eq.(^^ and the appropriate value of the 
cooling function Acooi is calculat ed. We based our numerica l 
code on the code described by Thoul fc Weinberg (1995), 



which is the standard, second-order accurate, Lagrangian 
finite-difference scheme. Details of the code are described in 



Stachniewicz fc Kutschera (2001) 



We handle timesteps and central bou ndary conditions in 
the way proposed by Thoul and Weinberg Thoul & Weinberg 



(1995). For dark matter, we treat the center as a hard sphere 
of some "small" radius re- In order not to affect the results 
of calculations the value of Vc should be much less than 
any other characteristic radius in the problem but it should 
not be too small because smaller means worse energy 
conservation and longer computation time. We have chosen 
Tc equal to the initial radius of the most innermost dark 
matter shell. For the baryonic component, if a shell falls 
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Figure 3. Test of accuracy - comparison of runs with various 
number of baryonic and dark matter shells 



below we assume that it has 'collapsed' and fix its radius, 
temperature and chemical composition. 

In the collapse of baryon matter one encounters the for- 
mation of shocks. In our code shocks are treated with the 



1967) 



artificial viscosity technique (R,ichtmyer & Morton 

It is worth to mention that set of equations |l6| is very 
stiff and may requir e very small timesteps . To deal with this 

and we used the 



problem we follow 



Haiman et al 



STIFBS routine from press et al 



(119961) 

mm. 



We have performed two sets of calculations, for the A- 
dominated Universe (Qi^h^ = 0.02, Q,m = 0.35, Qa = 0.65) 
and for the pure CDM model {Qth^ = 0.02, Qm = 1.0), with 
h — 0.72 in both cases. We started our calculations at z = 
500. The initial gas temperature and chemical composition 
were obtained by running our own program that starts at 
z = 10000 with equilibrium values and ends at a required 
redshift. We h a ve co mpared its predictions to the results of 
Galli & Palla (199J) and the differences were not greater 
than 10%. One should note note that these authors have 
included more species (e.g. deuterium and lithium) and some 
of their reaction rates were slightly different. The difference 
between the temperatures of the baryonic gas and the CBR 
turned out to be lower than 1% a,t z = 500. 

For the A-dominated model we have performed runs 
with cloud masses Mo from 2000 to 10^ M©, for the pure 
CDM model cloud masses were between 500 and 5000 M0. 
For each cloud mass we have chosen a few values of the initial 
dark matter overdensity between 0.07 and 0.40, depending 
on the mass. 

We have divided the baryonic component to 200 concen- 
tric shells of equal mass, the dark matter component was di- 
vided to 1000 equal mass shells. It turned out to be enough - 
in Fig. ^ we compare the calculated amount of virialized and 
collapsed mass for the same initial conditions and three dif- 
ferent runs: with 200 baryonic and 1000 dark matter shells, 
200 baryonic and 2000 dark matter shells and, finally, with 
400 baryonic and 2000 dark matter shells. The results are 
the same within our numerical accuracy. 
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6 RESULTS 

We expect that there wiU be a significant difference between 
the low-mass structure formation in the ACDM and the 
Standard CDM models. The main reason is that the rms 
fluctuations in the dark energy dominated model are about 
2.5-2.7 times lower than in the Standard CDM model (table 
|l|). The other one is that in the Standard CDM the average 
matter density is equal to the critical density so any overden- 
sity would eventually stop its expansion and collapse. In the 
ACDM the matter density is lower, so there is some thresh- 
old overdensity necessary for recollapse. Any given overden- 
sity would collapse slower (if at all) than in Standard CDM. 
It means that in the A-dominated Universe there would be 
much less luminous objects at large redshifts and the first 
ones would appear much later as compared to the dark mat- 
ter dominated Universe. In addition, we may expect that the 
masses of the first collapsed objects would be greater. 

The behaviour of dark matter shells plays the crucial 
role in the cloud collapse as there is much more dark mat- 
ter than the baryonic matter. By assumption, dark matter 
is coUisionless, so there is no energy loss mechanism apart 
from possible gravitational energy exchange. The evolution 
of the dark matter shells in the absence of baryonic matter 
is that their radii increase to some maximal value, then the 
shells collapse and after some oscillations due to bounces off 
the artificial hard sphere the dark matter cloud becomes sta- 
tionary. In the presence of baryonic matter the behaviour of 
dark matter shells is similar because the amount of baryonic 
matter is much lower than the amount of dark matter. The 
difference is that the average radius of a dark matter shell 
slowly decreases. 

The behaviour of baryon matter shells is that after 
reaching maximum radii, the shells collapse and the shock 
develops at so me stage due to press ure crowding of neigh- 
bouring shells (Haiman et al. 1996). The shock stops the 



collapse and the baryon shells virialize. Then if the clouds 
have a sufficient temperature and density for efficient cool- 
ing, they can undergo the final collapse. 

We assume that the shell has 'virialized' if the mean 
density inside the shell exceeds IStt^ times the actual mean 
baryon density of the Universe. This value was calculated for 
the top-h at overdensities in the Ein stein-de Sitter model but, 
however, Bryan & Norman (1998) have shown that in the 



ACDM models the dependence on Q,m is not very strong. Of 
course, this criterion depends on the redshift. For a virialized 
baryonic shell its radius and density are roughly constant 
but the mean matter density of the Universe decreases so 
the overdensity increases. This means that overdensity is not 
a good criterion of the collapse - we need to take something 
related to the actual size of the object. We have decided to 
assume that a shell has 'collapsed' if its radius falls below 
some 'small' value, in our case: the initial radius of the cloud 
-Ro. 

The results of our calculations are displayed in Figs. ^ 
- ^. Perhaps the most interesting are plots that show the 
evolution of the amount of the collapsed baryonic mass for 
various cloud masses and initial overdensities. Figs. ^ a-d 
show the fraction of the cloud mass that meets the collapse 
criterion at a given redshift for the A-dominated Universe. 
Similar plots for the dark matter dominated Universe are 
shown in Figs. |^ a-d. 
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Figure 4. The collapsed mass fraction as a function of redshift 
for Mo = 2x10^Mq, 3 x lO^M©, 2 x IO^Mq and IxIO'^Mq in the 
A-dominated Universe. Indicated are overdinsities corresponding 
to different curves. 
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If there is a line that is described but cannot be found 
on a plot it means that for this overdensity no baryonic shell 
has collapsed before z = 8. 

The results confirmed our expectations. In the dark en- 
ergy dominated Universe initial clouds need greater mass 
and larger initial overdensities in order to at least partially 
collapse before 2: = 8. If we restrict ourselves to clouds of 
z — 500 DM overdensities not greater than 0.4, the mini- 
mal mass needed for at least partial collapse before z = 8 
is 3000 M0 (Fig. 1^ b). For PI = 2OOOM0 and 5dm, i = 0.4 
the fraction of the collapsed mass is on order two percent 
(Fig. ^ a) . One should remember that for these cloud masses 
the rms ffuctuations at z = 500 are about 0.05 so overdensi- 
ties equal to 0.4 represent 8a peaks, which are extremely 
unlikely. For the dark matter dominated Universe, the min- 
imal cloud mass is about 500 M0 (Fig. ^ a) - for the initial 
overdensity equal to 0.4 (the rms overdensity is about 0.14) 
about 9% of the cloud mass may collapse before z = 8. 

If we look at higher masses, the collapse becomes easier 
and easier, threshold initial overdensities for partial collapse 
before 2 = 8 fall much faster than the rms overdensities. 
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Figure 6. Approximate redshift when the object starts its col- 
lapse. 
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Figure 7. Evolution of the baryonic (upper plots) and dark mat- 
tor (bottom plots) shells for Mo = 5 X IO^Mq, <5dm,i = 0.30 
for the A-dominated Universe (left) and for Mq = 5 x 10"^ Mq, 
^dm,i = 0.20 for the CDM-dominated Universe (right). 



However, for the A-dominated Universe and cloud masses 
not greater than 10^ M© the clouds need overdensities larger 
than 3(7 in order to at least partially collapse - unlike in the 
flat CDM model. For greater baryonic masses the collapse is 
a bit more likely but even for 10^ Mq (Fig. ^ d) they need 
at least 2 ct to start their collapse before z = 8. 

Fig. ^ shows approximately the redshift when the most 
innermost shell has collapsed for various masses and 3cr over- 
densities in the ACDM Universe and for 2a overdensities in 
the flat CDM Universe. 

In Figs. 1^ a-d we compare the evolution of dark mat- 
ter shells and baryon matter shells of the cloud with mass 
Mo = 5OOOM0, for the dark energy and dark matter dom- 
inated Universes. The values of initial dark matter density 
enhancements are, respectively, 5dm,i = 0.30 and 5dm.,i = 
0.20. 

The shells shown on the plots enclose 7%, 
17%, 27%, ...,97% of the bound mass. This fractional 
division of mass applies to both dark matter and baryon 
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Figure 8. Evolution of the shell temperatures and chemical evo- 
lution for a shell enclosing 12% of the bound mass, for Mo = 
5 X 10'^ Mq, 5dm,i = 0.30 for the A-dominated Universe (left) 
and for Mq = 5 X 10^ Mq, 5am,i = 0.20 for the CDM-dominated 
Universe (right). 



matter shells. It applies also to the Fig. |8| (plots a and b) 
that show thermal evolution of gas shells. 

The behaviour of both baryonic and dar k components 



i s quit e typical, very similar to the results of Haunanetal 
(ll996|) and our previous work ( Stachniewicz fc Kutschere 
pool] ). Dark matter shells expand to some maximum radius, 
then recoUapse to about a half of this value and then very 
slowly contract. In our calculations the behaviour of dark 
matter is quite independent of the cloud mass - however, 
the position and the value of the maximal radius depend on 
the initial overdensity. In contrast, the behaviour of baryon 
matter shells strongly depends on the cloud mass and initial 
overdensity. Of course, higher masses and overdensities cause 
faster and more violent final collapse. 

In Figs. ^ a and b we show the temperature of shells as 
a function of redshift. The curves indicating temperatures 
of baryonic shells behave in an opposite way than the radii 
of baryonic shells. The temperature falls during expansion 
and increases in the collapse phase. After virialization the 
temperature remains roughly constant until it starts to fall 
quite rapidly during the final collapse. These figures show 
that virial temperatures in Standard CDM are higher than 
in ACDM for objects of the same mass: even though for 
the ACDM cloud initial overdensity was higher (0.3 com- 
pared with 0.2), virial temperatures of its shells were about 
two times smaller. We also have noticed that in both cos- 
mological models shells that undego the final collapse must 
reach some critical temperature. For ACDM it is about 500 
K while for Standard CDM it is about 450 K and it is quite 
independent on the mass of the cloud. Shells at lower tem- 
peratures do not collapse before z = 8. 

The abundance (by mass) of various species as a func- 
tion of redshift is shown in Figs. |^ c and d. The results 
correspond to the shell of mass M = 0.12Mz of the same 
cloud as in Figs. |^ and |^ a, b. One can notice the increase of 
the amount of H2 molecules at later redshift when the shell 
collapses. This makes H2 cooling more efficient, triggering 
the final collapse of inner shells. 

In Figs. I^a-f we show plots similar to Fig. Qbut normal- 
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Figure 9. Normalized evolution of baryonic shells for A/q = 5 x 
^dm,i — 0.30 for the A- dominated Universe (left) and for 
Mo = SxlO^M©, 

^dmA — 0.25 for the CDM-dominated Universe 
(right) with cooling (upper plots), without H2 cooling (plots in 
the middle) and without any cooling or heating (bottom plots). 



ized to the maximum radius Rq and time of the turnaround 
to for Mo = 5 X IO^Mq, Sdmj = 0.30 for the A-dominated 
Universe (plots a, c and e) and for Mo = 5 x 10'^ M0, 
Sdm,i = 0.25 for the CDM-dominated Universe (plots b, d 
and f) with cooling, without H2 cooling and without any 
cooUng or heating (including photoioriization) . While plots 
with no cooling/heating (e and f) and no H2 cooling (c and 
d) are almost self-similar for different shells (the difference 
in final radii is probably due to different mean initial over- 
density for e ach shell) and sim ilar to those predicted by 
Bertschinger ( Bertschingerl 1985 ) , plots with cooling are self- 



nal collapse. In the A-dominated Universe perhaps low-mass 
objects may form through fragmentation of greater objects. 
Direct formation of such objects before z ~ 10 seems un- 
likely. 



REFERENCES 

Bardeen J.M. et al., 1986, ApJ 304, 15 
Barkana R., Loeb, A., 2001, ARA&A 39, 19 
Bertschinger E., 1985, ApJS, 58, 39 
Bryan G.L., Norman M., 1998, ApJ, 495, 80 
Eisenstein D.J., Hu W., 1999, ApJ 511, 5 
Fan X. et al., 2001, AJ 122, 2833 
Gain D., Palla P., 1998, A&A, 335, 403 
Haiman Z., Thoul A.A., Loeb A., 1996, ApJ 464, 523 
Katz N., Weinberg D.H., Hernquist L., 1996, ApJS, 105, 
19 

Lee A.T. et al., 2001, ApJ 561, LI 

Miller A.D. et al., 1999, A&AS 195, 5508 

Netterfield C.B. et al., 2002, ApJ 571, 604 

Padin S. et al., 2001, ApJ 549, LI 

Perlmutter S. et al., 1999, ApJ 517, 565 

Press W.H., Teukolsky S.A., Vetterling, W.T. Flannery 
B.P., 1996, Numerical Recipes in C. The Art of Scien- 
tific Computing. Second Edition. Cambridge University 
Press (1996) 

Pryke C. et al., 2002, ApJ 568, 46 

Richtmyer R., Morton K.W., 1967, Difference methods for 
Initial- Value Problems, Interscience Publishers, New York 
Riess A.G. et al., 1998, AJ 116, 1009 
Seljak U., Zaldarriaga M., 1996, ApJ 469, 437 
Stachniewicz S., Kutschera M., 2001, Acta Phys. Pol. B 32, 
227 

Thoul A. A., Weinberg D.H., 1995, ApJ 442, 480 



similar and Bertschinger-like up to virialization only. It also 
shows that H2 cooling is necessary to collapse before z — 8. 



7 CONCLUSIONS 

For the dark energy dominated Universe we need very large 
overdensities, usually much exceeding the 3cr limit, for low- 
mass (A/o < 10^ A/0) objects in order to at least partially 
collapse before z = 8. For the flat, dark matter dominated 
Universe the collapse is much faster and direct formation 
of a low-mass object is quite likely, even as soon as for 
2: ~ 20 — 30. It is a result of much lower normalization in 
ACDM (Table 1) and lower virial temperatures for ACDM 
clouds of the same mass and initial overdensity. A secondary 
reason may be that in the flat CDM model any overdensity 
will stop its evolution and possibly cool and collapse while 
for < 1 there is some threshold overdensity that de- 
pends on the model and epoch. The chemical evolution is 
pretty much insensitive to cosmology, at least before the fi- 



